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The paper presents some solutions to the five dimensional Einstein equations due to a perfect 
fluid on the brane with pure dust filling the entire bulk in one case and a cosmological constant (or 
vacuum) in the bulk for the second case. In the first case, there is a linear relationship between 
isotropic pressure, energy density and the brane tension, while in the second case, the perfect fluid 
\q . is assumed to be in the form of chaplygin gas. Cosmological solutions are found both for brane and 

■ bulk scenarios and some interesting features are obtained for the chaplygin gas on the brane which 

' are distinctly different from the standard cosmology in four dimensions. 
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I. INTRODUCTION 



Recently there has been a proliferation of brane world models, where the standard model matter is confined 
in a four dimensional space time said to be a singular hypersurface or 3-brane embedded in a (4 + d)-dimensional 
space time. This (4 + d)-dimensional space time is said to be the bulk. Though the matter is all confined 
to a 3-brane, the gravitational field can propagate in the higher dimensional bulk. The d extra dimensions, 
however, need not be small or even compact in this theory. In fact, Randall and Sundrum [1] have shown 
that for d — 1 the fifth dimension may be infinite. The effective equations for gravity in four dimensions were 
obtained by Shiromizu et al [2] using Israel's boundary conditions and ^-symmetry for the bulk space time, in 
which the brane is embedded. It has been demonstrated by Binetruy et al [3,4] that generically the equations 
governing the cosmological evolution of the brane will be different from those corresponding to the analogous 
Friedmann equations of standard cosmology. Essentially the difference lies in the fact that the energy den- 
sity of the brane appears in quadratic form in the brane equations, whereas it is linear in the standard cosmology. 

In the present paper our motivation is to solve the five dimensional Einstein equations for a perfect fluid 
• • . on the brane while the matter in the bulk is cither in the form of pure dust or a cosmological constant (or 
vacuum). It is well known that in the standard cosmology, a special kind of fluid called chaplygin gas [5 
and references therein] has a very interesting property that it can drive the standard AD universe from its 
initial dust dominated stage to a final A dominated stage at the end, so that at the late stage, one can obtain 
an accelerating world model consistent with the recent Supernova type la observations [6]. Measurements 
of luminosity-redshift relation, obtained from observations suggest that most of the energy in the universe 
happens to be dark. Different kinds of the dark energy viewed from different angles have so far been proposed 
[7] to explain the late time acceleration of the universe. Recently, Ponce de Leon [8-10] studied the acceleration 
of the 4D universe as a consequence of the time evolution of the vacuum energy or the brane tension in the 
context of the brane world models. 



In our paper, section II presents the derivation of brane equations starting from the Einstein equations 
in the bulk with a general form of matter there. Then in section III, assuming pure dust in the bulk, 
Einstein equations in the brane show that only a kind of barotropic equation of state pb — — |pj holds 
between the effective isotropic pressure pb and effective energy density pb (i.e, a linear relation among the 
actual pressure and energy density in the brane and the brane tension). Cosmological solutions are obtained 
for both bulk and brane scenario. Section IV deals with perfect fluid in the form of chaplygin gas in the 
brane model. When the bulk is vacuum, then the dynamics shows initial deceleration of the brane with 
subsequent acceleration going over to the ACDM model. On the other hand, when the bulk has negative 5D 
cosmological constant, i.e, an anti de-Sitter bulk, there occur two different cases. Depending on the param- 
eters in the field equations, the brane may either transit from a decelerating stage to an accelerating stage 
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in course of time or it may simply exhibit a recollapsing behaviour. The paper ends with the concluding remarks. 



II. BASIC EQUATIONS 

We shall assume that the geometry of the five dimensional bulk is characterized by the space time metric of 
the form 

ds 2 = g^dx^dx" = g l:j dx l dx j + b 2 dy 2 (1) 

with y as the fifth coordinate. The hypersurface y = is identified as the world volume of the brane that 
forms our universe. As usual, for simplicity, the explicit form of the metric (1) is taken to be 

ds 2 = -n 2 (t,y)dt 2 +a 2 (t,y)5 l:j dx l dx : > +b 2 {t,y)dy 2 (2) 
where for simplicity we choose the usual spatial section of the brane to be flat. 

Now for cosmological solutions, we shall have to specify the matter explicitly. In fact, we shall classify the 
matter as of two distinct parts namely 

i) matter confined to our brane universe and 
ii) matter distributed over the 5D bulk. 

Hence, the energy momentum tensor can be decomposed into 

rpV rpV I rp V ( c \\ 

where the explicit form of the matter is [11, 3, 4] 

f% bulk = diag(-p B ,p B ,PB,PB,P5) (4) 

and 

V brane = — diag(- Pb , Pb , Pb , Pb , 0) (5) 

In general, energy-momentum tensor in the bulk depends on t and y (i.e, pb = PB(t,y)) but to recover a 
homogeneous cosmology in the brane, matter field in the brane is assumed to be the function of time alone. 
Also the brane is assumed to be infinitely thin. Now the explicit form of the five dimensional Einstein equations 



for the above metric (2) having matter distribution in the form (3) are [11, 3] 
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where an overdot stands for differentiation with respect to time £ and a prime indicates derivative with 
respect to the fifth coordinate y. 

For brane cosmology, we need solution of the Einstein's equations (7)-(10) as y — > 0. Now, for a well defined 
geometry, the metric tensor should be continuous across the brane (i.c, y — 0) while its derivative with respect 
to y may be discontinuous at the hypersurface y = 0. As a result, Dirac delta function will appear in the 
second derivatives of the metric coefficients with respect to y namely [3] 



a" = (a") + [a'] 6(y) 



(11) 



Here (a") is the non-distributional part of the second derivative i.e, the usual second order derivative and [a'] 
is the measure of discontinuity of first derivative across y — (with S(y), the Dirac delta function) as 



[a'} = o'(0+) - o'(O-) 



(12) 



If we now proceed to the limit as y — > and match the Dirac delta function in both sides of the field equations 
(7) and (8) we obtain 



ao&o 



(13) 
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y (3p b + 2p b ) 



(14) 



where subscript '0' in the scale factors indicate their values on the brane. 



One may note that the jump of the field equation (9) using equations (13) and (14) results in the energy 
conservation relation on the brane namely 



pb + 3 (pb+Pb) — = 
a 



(15) 



Further, if we take the average value of the field equation (10) for y — > + and y — > and impose the 
Z 2 -symmetry we get [3] 



d'o do 2 k 4 , . „ x « 2 P5 

1" — 2" = PHPfc + 3 P&) 



a 
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(16) 



Here we have chosen no = 1, which is allowed by a suitable time transformation and hence t is now the 
usual cosmic time on the brane. This equation can be termed as generalized Friedmann type equation on 
the brane. The basic difference with the corresponding equation in standard cosmology is that, here the square 
of the Hubble parameter depends quadratically on the brane energy density, in contrast with the usual linear 
dependence for Friedmann universe. 
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Moreover, by introducing a function £(i, y) denned by 

«'.»)- ("> 

the five dimensional field equations (7) and (10) can be written in compact form as (for y ^ 0) 

= -| a'a 3 K 2 p B (18) 

£ - | da 3 K 2 p 5 (19) 
Now equating the time derivative of equation (18) with the y-derivative of equation (19) we have 

a' P B + ap' 5 + (p B + p 5 ) (V + 3^ = (20) 

As we have assumed T05 = indicating that there is no flow of energy along the extra dimension, the Bianchi 
identity V^G^° = results 

Pb + ^(p b +Pb) + ^(pb+P 5 ) =0 (21) 

Also taking average value of the equation(17) for y — > + and y — > 0~, imposing the .^-symmetry and using 
the junction condition (13) we obtain the generalized Friedmann equation 

where ^0 is the value of £ on the brane. 

Therefore, the cosmological evolution on the bulk is completely characterized by the Einstein field equations 
(7)-(10) (also (18), (19)) and by the conservation equations (20) and (21). On the other hand, brane cosmology 
will be determined by equations (16) and (22), junction conditions (13) and (14) and by the matter conservation 
on the brane (15). 



III. BRANE COSMOLOGY FOR DUST FILLED BULK 

If the bulk matter is assumed to be in the form of pure dust i.e, pb — P5 = 0, then integrating the conservation 
equation (21) we get 



PB(t,y) = Pl (y)/a 3 b (23) 

where pi(y) is function of the coordinate y alone. 
Also comparing the conservation equations (20) and (21) and using the field equation (9) we have 



n' = (so that n may be assumed to be unity without the loss of generality) and a' = a(y) b (24) 
with a, an arbitrary function of y. 
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As P5 = 0, the equation (19) yields £ = /3(y), a function of y alone which again in view of (18) satisfies the 
relation 

fty) = -!?f£tfy (25) 

Now using (24) in (17), we obtain the differential equation for a in the form 

a 2 a 2 = a 2 a 2 — j3 

which on integration gives the solution for the scale factor on branc in the form 

a 2 (i,y) = -i-j [fi{y) + {a 2 (y)t + t (y)} 2 ] (26) 
with to(y), an arbitrary integration function. 

In order to obtain cosmological solutions on the brane, we now consider the generalized Fricdmann equation 
(22) written earlier 



do 2 _ 2 _ 00 
al ~ 36 Pb a; 







Further, combining (24) with the junction condition (14) we get for the effective equation of state for the brane 
matter 



Pb = -g Pb (28) 



Equivalently, as 



p b = p - A and p b = p + A (29) 
with A representing the brane tension, the above equation of state yields the following linear relation 

P=-\p+\>> (30) 

Thus for dust bulk, the fluid in the brane must satisfy the above linear relation between isotropic pressure, 
energy density and brane tension. Using now the above equation of state (28) into the brane conservation 
equation (15) we get 

Pb = Po/ a o j Po here is simply a constant. (31) 
Now since a' — a(y)b as is evident from the equation (24), the brane must satisfy the following relation 

[a'] = 2b a 

and hence using equations (13) and (31), the arbitrary constants po and an above are shown to be related as 

a = -K 2 p Q /6 (32) 

Further, for the dust model in the bulk and the equation of state (28) in the brane, the Friedmann type 
equation (16) simplifies to 



G 



K 4 2 

2aoa'o + a 'o 2 = M > where /j, = (33) 

18 



which has the first integral 



4>' = S + 7* (34) 



• 2 _ M , Mo 

where is an integration constant. This is fully consistent with the equation (27) provided /Uq = — 0o 



Now further integration of the equation (34) yields the explicit expression for the brane scale factor do as 

a 2 a = 2/i (ti ± tf - (35) 
A 4 

Here since d'o = — we always get either an accelerating model or a decelerating model of the brane 
universe depending on the sign of the constant fj,o, but there can not be any midway transition from the 
decelerating phase to the accelerating one (or the reverse). The choice of sign in front of t in equation (35) 
leads to an expanding (+vc sign) or a contracting (-ve sign) brane world. 



IV. CHAPLYGIN GAS BRANE MODEL 

We consider here a modified chaplygin gas [12], which satisfies an equation of state in the following form 

Pb =Ap h -\ (36) 
Pb a 

where A and B{> 0) are constants and < a < 1. 

So in this case from the energy conservation equation (15), the expression for matter density is given by [12] 



Pb 



1 > 92 +B 



(37) 



Ui I 3(l+A)(l+a) 

with p 2 , an arbitrary integration constant. 
Case-I: Vacuum Bulk 

As there is no matter in the bulk we have T" bulk = and the field equations (18) and (19) show that £ must 
be a constant (say, C). Hence we have the generalized Friedmann equation (see eq. (22)) 



Also the generalized Friedmann type equation (16) now simplifies to 

- + i = ~Pb(Pb + 3p b ) (39) 

Further, our calculation will be much simplified if we choose a = 1 in the equation of state (36) and hence 
from (37) we get the differential equation for the scale factor a in the form 
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a a = - 



k 4 (3A + 2) po 
36 (A+l) a (6A +5 ) 



+ 



36(1 + A) 



(40) 



/ \ 1/2 

Note that as a — > (early epoch) p fc — > oo, d' — > — oo while as a ^ oo, p b — > ( ^4 J which in view of 



-1/2 

(36) leads to pb — > — ( J and d' — ► 00. This shows that asymptotically the equation of state reduces to 

Pb = —pb, which points to a ACDM model in the corresponding FRW universe of standard cosmology. But 
for some special cases of the brane world models (shown in what follows) the situation is different, since here 
d'o may remain negative throughout even if the scale factor ao increases indefinitely Such behaviour is not 
consistent in the four dimensional Friedmann model containing chaplygin gas, which always asymptotically 
goes over to the ACDM model with the increase of the scale factor. 



We now proceed to solve equation (40) for A = | (for which, solution in closed form is possible). The first 
integral of equation (40) can be written in an integral form as 



dz 



\Jbz 2 + Cz + d 



= ± (t-to) 



(41) 



with z = <2q, b = d = -^jf 2 -, solutions of which are given by 



a n = 



VAbd - C 2 
~2b 



Sinh 



±4Vb (t - t ) 



(when Abd - C 2 > 0) 



(42) 



and 



VC 2 - Abd 
2b 



Cosh 



±4Vb (t - t ) 



4 = Yb 
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(when 46d - C 2 < 0) 



(when Abd - C 2 = 0) 



(43) 



(44) 



Graphically, it can be seen that for equation (42), the scale factor expands exponentially whereas for equation 
(43), it decreases to a minimum before further expansion. 

Case-II: Bulk with Cosmological Constant 

In this case we choose in the bulk pb = Pb = —pB- It follows then from Bianchi identity, i.e, the conservation 
relation (21) demands ps = 0, which in turn in view of (20) yields p' B = 0. Hence ps is purely a constant said 
to be the cosmological constant in five dimension (= A5). Now the generalized Friedmann like equations (38) 
and (40) are modified in this case as 
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(46) 







We choose here A — 1/3, so that at the begining, the brane world may be said to be filled up with (p = p/3) 
radiation. 
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Fig. 1(a) 



Fig.l(b) 



Fig. 1(a) shows the variation of the scale factor (given by eq.(42)) over the time for k = 2, B = po = 3, C — 1 and 
to = 1 while Fig. 1(b) shows the same (for eq.(43)) for n = 2, B = po = 3, C = 3 and to = 1 In both the cases, (+)ve 

sign is taken. 



We now examine the following cases: 



For this restriction d'o < for all time while o'o vanishes at some finite do- It is therefore, a decelerating 
model of the universe with only a maximum but no minimum. It is a recollapsing brane model (see fig. 2(a)) 
with ao — > and d'o — > — oo at both ends (radiation to radiation). 



Here Hq is non-zero throughout the evolution. As ao — ► 0, d'o — > — oo and as ao — > oo, d'o — ► oo. So there 
is a transition from deceleration to acceleration in course of evolution of the brane (Graphically it is shown in 



In this case there is always d'o < with Hq — > as ao — > oo. So this has a typical Einstein-de Sitter like 
behaviour for the spatially flat standard model (presented in fig. 2(c)). The interesting feature of such a brane 
world model is that even though at large a the fluid has an equation of state pb = —pb, the model does not 
show the character of a ACDM model as in standard cosmology. Rather it still remains a decelerating model. 
As we note that in one of the above three models of the brane universe, for example in the second model, 
there occurs a transition from a decelerating to an accelerating phase. It may be worthwhile to locate the said 
transition point where d'o = 0. It is, however, not difficult to find out from the equation (46) that the brane 
world model with A = 1/3 remains accelerating (d'o > 0) as long as 



a) When ^-^|A 5 | <0 



b) When^-^|A 5 |>0 



fig. 2(b)). 



c) When^-^|A 5 |=0 i.e, |A 5 | = ^ 




that is 




which finally reduces to the following condition by a straightforward calculation: 




In the standard cosmology with the chaplygin gas in FRW universe, the corresponding inequality is 
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50 100 150 200 Fig. 2(c) 



Fig. 2(a) shows the recollapsing model of the universe for k = 1, B — 1, p = 48, C = 1, |As| = 1 and t = 24. 
Fig. 2(b) shows a transition from deceleration to acceleration for k = 1, B — 16, po = 24, C = |, |As| = 1 and to = 1- 
Lastly, Fig. 2(c) shows the Einstein-de Sitter universe for k = 1, B — 8, po = 48, C = 1, |As| = 1 and to = 1- 



a s > po/B. 

Note that, this is only a formal comparison, as p , B etc. in the brane and the corresponding quantities in 
FRW universe are not really identical. 



V. CONCLUDING REMARKS 

When the bulk is filled with only dust, there are certain restrictions imposed on the metric as well as on the 
matter content of the brane. As results of such restrictions the metric coefficient g u becomes independent of 
the y-coordinate, while the perfect fluid in the brane satisfies an effective equation of state pb = —^Pb, which 
apparently indicates that matter in the brane is in the form of dark energy. 

On the other hand, if the matter in the bulk is in the form of cosmological constant or vacuum, there is 
neither any restriction on the metric nor on the matter contained by the brane. Binctruy ct al and Leon 
[3, 4] have shown perfect fluid solution with barotropic equation of state for brane matter. In this paper, 
we have shown that for brane matter, it is possible to have chaplygin gas form of fluid and obtained some 
interesting conclusions. The recollapsing model of the brane universe shows that in course of evolution it 
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reaches a maximum size and finally collapses to a singularity. So for such a model, there is deceleration 
through the entire evolution period and it has singularities at both ends. This difference in behaviour of the 
model is caused by the modifications introduced in the field equations by the presence of an anti-de Sitter bulk 
(A5 < 0). Another point of significance for a chaplygin gas in the brane world model is that, the density does 
not vanish as the scale factor increases indefinitely, rather it reaches a finite magnitude with an equation of 
state pb = —pb- This behaviour still remains valid in one of the brane world models, which expands indefinitely 
with deceleration throughout and resembles the standard Einstein-de Sitter model. Such behaviour of the brane 
is in clear contradiction to the standard AD Fricdmann universe. 
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